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Abstract. We present a discussion of the charge response and the charge 
fluctuations of mesoscopic chaotic cavities in terms of a generalized Wigner- 
Smith matrix. The Wigner-Smith matrix is well known in investigations of time- 
delay of quantum scattering. It is expressed in terms of the scattering matrix 
and its derivatives with energy. We consider a similar matrix but instead of 
an energy derivative we investigate the derivative with regard to the electric 
potential. The resulting matrix is then the operator of charge. If this charge 
operator is combined with a self-consistent treatment of Coulomb interaction, the 
charge operator determines the capacitance of the system, the non-dissipative 
ac-linear response, the RC-time with a novel charge relaxation resistance, and 
in the presence of transport a resistance that governs the displacement currents 
induced into a nearby conductor. In particular these capacitances and resistances 
determine the relaxation rate and dephasing rate of a nearby qubit (a double 
quantum dot). We discuss the role of screening of mesoscopic chaotic detectors. 
Coulomb interaction effects in quantum pumping and in photon assisted electron- 
hole shot noise are treated similarly. For the latter we present novel results for 
chaotic cavities with non-ideal leads. 
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1. Introduction 

Quantum transport in structures so small that the quantum wave nature of particles 
becomes important presents many theoretical and experimental challenges. The 
physical realization of such small structures is quite diverse and ranges from electronic 
transport through (partially) coherent samples [1, 2, 3, 4, 5] to scattering of photons 
in cavities [6, 7] or scattering of particles in compound nuclei [8]. Often the dynamics 
of particles inside the scattering region is chaotic. In an optical cavity this is a 
consequence of the irregular shape of the resonator and, similarly, in electrical samples 
chaotic scattering results from the boundary or the impurity configuration. This 
paper, devoted to electron transport, is focused mostly on the particular example of 
a chaotic cavity, namely a two-dimensional chaotic quantum dot. A typical quantum 
dot is usually patterned into the two-dimensional electron gas at the interface of semi- 
conducting heterostructures. The two-dimensional electron gas results from strong 
quantization of the electron motion perpendicular to the interface. The shape of the 
sample is formed with the help of metallic gates on top of the structure. Even if the 
dot is free of impurities, the shape of the dot is usually quite irregular, so that the 
dynamics of electrons is classically chaotic. (If in addition there are impurities chaos 
results independently of the geometric shape of the dot). 

At sufficiently low temperatures, when the length over which the carriers retain 
phase memory becomes comparable to the dimensions of the structure we enter the 
regime of mesoscopic physics [9, 10]. In this regime, every dot exhibits fluctuating 
physical properties due to the high sensitivity of quantum interference effects on the 
exact geometry of the dot. Then properties of small structures must be described by 
their statistics in the form of mesoscopic (sample-to-sample) distribution functions 
rather then by their averages over ensemble of similarly fabricated samples. 

We are interested in dots which are connected to the outside via one or several 
leads which permit the exchange of carriers with electrical contacts (reservoirs). The 
connection between the dot and reservoir can be highly transparent or alternatively 
we can insert tunnel barriers which increasingly insulate the dot from the contacts. 
If the dot is closed, its equilibrium charge changes in response to the voltage applied 
to the external gates. Except for special gate voltages the charge on the dot does 
not fluctuate. For poorly transmitting contacts (strong tunnel barriers) the charge of 
the dot is strongly quantized and at low temperatures this quantization can block 
transport (Coulomb blockade). Such a blockade is important if a typical energy 
needed to add an electron into a dot, ~ e 2 /C (capacitance C defined by geometry), 
is large compared to the temperature and the escape rate 7 OS c of carriers from the 
dot, e 2 /C ^> fceT, Hj esc [9]. However, when the barrier becomes transparent and the 
contacts are wide open the charge quantization vanishes and weak charge fluctuations 
appear at all voltages. 

The main subject of this paper is the role of Coulomb interactions on transport 
properties of the dots which are connected with highly transparent contacts to 
reservoirs. It turns out, that in many transport problems the role of Coulomb 
interactions is closely connected to the Wigner-Smith matrix [11, 12]. This matrix 
is well known in investigations of the time-delay of quantum scattering [13, 14, 15, 16, 
17, 18]. It is expressed in terms of the scattering matrix and its derivatives with energy. 
We consider a similar matrix but instead of an energy derivative we investigate the 
derivative with regard to the electric potential [19, 20]. The resulting matrix is then 
the operator of charge. If this expression for the charge operator is combined with 
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a description of interaction on the Hartree level (random phase approximation) [21] 
the charge operator thus found plays an important role in a number of transport 
problems. The list of these problems includes charge rearrangement in weakly 
non-linear transport, the current-response to oscillating potentials, relaxation and 
dephasing in weakly coupled nearby conductors, the theory of quantum detectors, 
adiabatic quantum pumping, and frequency-dependent thermal and shot noise. 

There are several important questions one might ask: What are the signatures 
of Coulomb interactions in transport properties of a sufficiently open quantum dot? 
Which experiment could probe the Coulomb interactions? Can Coulomb effects be 
distinguished from other effects like dephasing? Of course for poor contacts the 
Coulomb blockade has a very clear signature in transport. However, many experiments 
are performed with few- and multi-channel open dots, so it is also important to account 
for Coulomb interactions in such samples. 

Our discussion is based on the scattering matrix approach to chaotic systems 
[22, 23, 24]. The formalism is a subject of many extensive reviews [2, 3, 5, 9], and we 
just outline here why this method is advantageous for our purposes. Our starting point 
is the energy-dependent scattering matrix S(e) that relates incoming and outgoing 
electrons at energy e. This matrix is obtained from the solution to the Schrodinger 
equation for non-interacting electrons. For low frequencies ui <C 7 CS c, the scattering 
matrix varies weakly on the scale ~ to, so for linear transport one can expand transport 
coefficients in ui. It is this expansion which makes the Wigner-Smith matrix appear. 
For quantum dots the statistical properties of the Wigner-Smith matrix have recently 
been widely explored [13, 14, 15, 16, 17, 18, 25], so that one can use them to find the 
statistical distribution of measured quantities. 

The scattering matrix approach is not the only possible method to solve transport 
problems. For a discussion of alternative methods and applicability of our approach 
we refer to Sec. 9. The self-consistent treatment of the Coulomb interaction is the 
subject of our paper, and we apply this method even for few-channel quantum dots. 

We proceed as follows: In Section 2 we first consider non-interacting electrons and 
introduce the Wigner-Smith matrix as a natural building block for a low-frequency 
low-temperature transport theory of open cavities. On the example of a single-channel 
cavity we demonstrate that taking into account the Coulomb interactions is essential 
in finite frequency transport. A self-consistent treatment allows us to consider the 
effect of Coulomb interactions and find the internal potential and its derivative with 
respect to external perturbations for a dot with arbitrary number of channels. Later 
we concentrate on applications of our approach to various experimental set-ups. Sec. 
3 considers the role of Coulomb interactions in the renormalization of the effective 
capacitance and introduces the mesoscopic charge relaxation resistance in terms of the 
Wigner-Smith matrix. These results are used later in Sections 4-8. Sec. 4 considers 
the displacement current induced into a nearby gate due to the charge fluctuations 
associated with shot-noise in a dc-biased two-channel cavity and finds its mesoscopic 
distribution. In Sec. 5 the relaxation rate and dephasing rate of a charge qubit near a 
chaotic cavity are discussed. In Sec. 6 we describe how to use these results to describe 
a quantum detector. Sec. 7 concerns 'quantum pumps', which pump electrons due a 
time-periodic variation of the dots shape. The effect of Coulomb interactions on the 
pumped voltage is considered self-consistently. As another application of the Hartree 
treatment we consider an ac-biased multi-channel quantum dot in the Sec. 8 and find 
the photon-assisted shot-noise both for completely open dots and partially open dots 
with equal channel transparencies. In this section we use results of Appendix A for 
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various correlators of the scattering matrices. We discuss other approaches and the 
applicability of our results in Sec. 9 and conclude in Sec. 10. 

2. Generalized Wigner-Smith matrix 

In this work we consider electrical transport problems which can be described with 
the help of a generalized Wigner-Smith matrix [11, 12]. We illustrate these transport 
problems for chaotic cavities. In linear response to oscillating potentials the Wigner- 
Smith matrix describes the leading non-dissipative response and in certain cases 
also the leading dissipative term. In particular, capacitances [26], charge relaxation 
resistances [26], and a Schottky resistance in the presence of shot noise [27], can be 
expressed in terms of interaction constants and the Wigner-Smith matrix. 

To be specific, consider the structure in Fig. 1. A cavity is connected with a single 
lead to a "contact". At the Fermi energy there are N open quantum channels. As a 
consequence there exists a scattering matrix S of dimension Nx N with elements S mn 
which relates the current amplitude of the incoming state n to the current amplitudes 
of the out-going state to. Each carrier incident on the cavity will after some delay 
be reflected and leave the cavity to be scattered back into the contact. According to 
Wigner and Smith the time carriers with energy E spend in the cavity can be found 
from the matrix, 

2iri dE 1 ' 

Taking the trace and multiply with Planck's constant gives a sojourn time t s = 
htvAf/N. The time defined in this way depends on how exactly we define the 
scattering matrix S. We could define S right at the entrance to the cavity (right most 
broken line in Fig. 1), somewhere in the long lead, at the entrance to the contact, or 
even deep inside the contact (left most broken line in Fig. 1). In the original work 
[12] the 5-matrix is clearly defined in an asymptotic way on the surface of a sphere 
with a radius that is eventually taken to infinity. The time calculated by Smith is a 
delay that is a difference of a time in the presence and absence of a scattering center. 
This is a sensible procedure if we are not interested exactly where a carrier incurs the 
delay. 

Our aim, however, is to obtain an expression for the charge, not the time-delay. 
Moreover, we are interested in the local charge. We have to know, whether this charge 
is accumulated inside the cavity, in the lead, or inside the contact. Therefore, the 
Wigner-Smith matrix is not the appropriate object [19, 20]. The local question which 
we are asking can of course also be asked for time. For instance, the widely discussed 
question, "How much time does a carrier spend while traversing a tunnel barrier?" 
is a question about local properties [19, 20] which can not be answered by an appeal 
to the Wigner-Smith matrix. It was already noticed in the treatment of this problem 
[28], that instead of using derivatives with regard to energy to find delay, we have 
to consider derivatives with regard to the electric potential U to find its conjugate, 
charge. Therefore, consider a small change of the potential inside the cavity only. We 
can then consider the matrix 

2m dyeil) 

Here the S-matrix can be defined asymptotically, i. e. in the contact to the left of Fig. 
1. Nevertheless, the answer given by Eq. (2) depends now only on the effect which a 
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Figure 1. Mcsoscopic capacitor connected via a single lead to an electron 
reservoir. V\ is the voltage applied to the contact, U is the potential in the 
cavity. 



small change of potential inside the cavity generates. We can now introduce the dwell 
time [20], which is the total integrated density in the cavity in an energy range dE 
divided by the incoming carrier flux (N/h)dE 

h , , h „ + dS 

T * = N trM = l N t ^d{euy ® 

We can choose any other potential, in the lead, or in the contact, and ask in a similar 
way what the dwell time is in the region of interest. We can take the region to be 
arbitrarily small and instead of a derivative of the 5-matrix investigate the functional 
derivative with regard to the local potential [29, 30, 31]. For the purpose of this work, 
it is, however, sufficient to consider derivatives with regard to single potential U. 

Instead of time-delay, or more precisely dwell time, we alternatively use the 
language of density of states. The density of states inside the cavity is given by 



v = txM=-^-M =-^-^ln(detS), (4) 

2m V d(eU) J 2m d(eU) y ' w 

where U is the potential in the cavity. If v is multiplied by an energy interval dE and 
the charge e, then AQ = evdE is the charge in the cavity due to the incident scattering 
states in the interval dE. Importantly this is the charge of non-interacting carriers. A 
conductor, when charged will respond with an induced electrical potential U to bring 
the charge in a given region to the value permitted by the Coulomb interaction. 

Suppose that we consider the linear response of the current in the contact 
dl(t) — dl{u) exp(-iujt) to an oscillating potential applied to the contact of the 
sample, dV(t) = dV(u) exp(-iujt). In the zero-temperature limit the response is to 
leading orders in frequency, 

G{uj) = dI{uj)/dV{uj) = -iuehiAf + uj 2 (e 2 /2h)tTAf 2 + 0{u 3 ). (5) 

The first term is a non-dissipative response determined by the density of states v 
of the cavity. e 2 u has the the dimension of capacitance and is called the quantum 
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Figure 2. Mesoscopic capacitor connected via a single lead to an electron 
reservoir and capacitively coupled to a gate. V\ and Vo are the potentials applied 
to the contacts, U is the electric potential of the cavity. 

capacitance [26, 32, 33, 34], C q — e 2 v. The second term is real and thus represents 
the leading dissipative part of the response. 

Equation (5) is not, however, a physically acceptable result already for the 
geometry shown in Fig. 1. It is necessary to consider the fact that electrons are 
interacting particles. To see the problem, suppose that the cavity and/or the lead are 
charged above what the ionic background of the conductor permits. As a consequence 
there will be long range electrical field lines emanating from the cavity and the lead. 
However, if the contact, the lead, and the cavity are the only metallic bodies, the 
electrical field lines emerging from the cavity or the lead will eventually end on the 
metallic contact. There exists, therefore, a Gauss volume taken large enough such 
that the electrical flux through its surface vanishes. The charge inside this Gauss 
volume is a constant of motion. Instead of Eq. (5) we should find G(w) — 0. We can 
see the profound effect of the long range Coulomb interaction in another way. The 
particle current density I p (r,t) is not necessarily a conserved quantity. But if we add 
the Maxwell displacement current 7 d (r,i) = eLdE(r,t)/dt with E(r,t) the electric 
field and consider the total current 7 tot (r, t) = I p (r, t) + I d (r, t) we have a quantity of 
vanishing divergence, div J tot (r,i) = 0. The total current has no sources or sinks. It 
is important to emphasize that it is the total current that is measured in experiment. 
Thus Eq. (5), valid for non- interacting carriers, is on physical grounds unacceptable. 

We now consider a more general arrangement and explain how the self-consistent 
Coulomb potential is determined. We consider transport through a dot with several 
contacts at which voltages V a (t) — V a cos(u;t + (f> a ) at some frequency ui are applied. 
In the end our results are also applicable for dc transport, if the limit lo — > is taken. 
In fact, in a dc-biased system, within linear transport, the potential U(t) = const 
and is not important, since at zero frequency, w = 0, the current is gauge invariant 
and knowledge of the internal potential is not necessary. However, it's derivative with 
respect to an external perturbation in the a-th lead, dU/dV a , can be used, e.g. in the 
analysis of the magnetic field asymmetry of the current in the non-linear transport 
through a dot [35]. As described above, the charge inflow into the dot shifts the 
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internal potential U(t) due to capacitive coupling to the gate C, kept at a potential 
V (t) = V cos(ut + <}> ). 

The current J a (fi) at finite frequency £1 in the a-th lead attached to the dot is 
the sum of two contributions: particle current i£(fi) and displacement current /^(fi). 
The former is due to variations of electro-chemical potentials in various leads and is 
expressed via conductances G aj g(fi) [26] of non-interacting electrons 

G af3 (n) = j Jdetr (5 af3 l a - l a 5+ (e)l p S(e + Ml)) /(g) ~ ^ + m) . (6) 

Here we introduce the matrices l a which have unit elements along the diagonal of 
the channels of contact a and zero otherwise. The displacement current, the result of 
screening, corresponds to variations of charge inside the dot and is expressed through 
the potential U(Q) and yet unknown susceptibility x Q (fi). Thus the current at probe 
a is 

la(fi) = I»{Sl) + = E G a p{H)Vp{H) + XaimW- (7) 

We require gauge-invariance, that implies that the global shift of all potentials 
Vi(fi) -> ^(Sl) - f(Q),U(n) -> ?7(0) - /(fi) by an arbitrary /(fi) docs not change 
the iotaZ current [26]. This determines the susceptibility x«(^) = — S/3 G a p(fl). As 
a consequence the current depends on difference of potentials only, as expected: 

4(0) = Ga/j(fi)W»(n) - W)- ( 8 ) 

On the other hand, from charge conservation ^ a ia(t) = C(d/di)[Vo(t) — U(t)], we 
find that Uq, ^ only at O = ±o>. In response to potentials with Fourier components 
y a ,Lo = V a exp(i(j) a ) the frequency-dependent potential U u is 

E Q /3 G "/3(w) +iwC 

At low frequencies, which is often important, the scattering matrices S are weakly 
dependent on energy, and the first terms in an expansion of the scattering matrix in 
energy are sufficient to evaluate the potential U w and the partial derivative of U u with 
respect to the potential V a ^. With the Wigner-Smith matrix Af introduced in Eq. 
(1) and the matrix of voltages V u = J2 a ^-ctV a ,u> the sample-specific expressions read 

tt _T/ _ trA/"(K, - \V .„) dU„ _ trJVl a 
u C/e^ + trM ' 8V a ^ ~ C/e 2 + trAA- ( ' 

Both quantities naturally vanish in the non-interacting limit C — > oo and reach a finite 
value in the opposite, strongly-interacting limit, C — > 0. In the following we use this 
approach to calculate several transport properties. 



3. Capacitance and charge relaxation resistance 

Fig. 2 shows a single-lead cavity which is separated by an insulating material from a 
back gate at voltage V (t) . Now we have the possibility to drive an ac-current through 
the system. In addition to the current response to V\ (t) we can investigate the response 
to Vo(t). A consistent treatment of the long range Coulomb interaction demands that 
these two responses are different only in sign. In fact the response must depend only 
on the voltage difference V(t) = V\(t) — Vo(t), see Eq. (10). Here we describe the 



Charge Flucutations 



8 



Coulomb interaction in simple terms using a geometrical capacitance C which links 
the potential difference between cavity and gate and the charge on the cavity and the 
gate. As a consequence the conductance of the cavity capacitively coupled to a gate 
is now given by 

G{lo) = -iwC^ + uj 2 C 2 R q + 0(u! 3 ). (11) 

Here the first term is the electrochemical capacitance [26] which is the scries 
capacitance of the geometrical capacitance C and the " quantum capacitance" ve 2 , 

C-^C^ + (ye 2 y\ (12) 

Thus the capacitance is not a purely geometrical quantity but depends via the density 
of states v on the properties of the specific electrical conductor. Notice, that C^ 1 is 
a fluctuating quantity, with mesoscopic average equal to (1/C„) = 1/C + A/e 2 . Here 
A is the mean level spacing in the dot. Equation (12) has important implications: for 
instance, if the cavity is deformed into a ring with an Aharonov-Bohm flux through 
the hole of the ring, the capacitance exhibits quantum oscillations which are periodic 
in the flux [36, 37]. Such oscillations have been observed in arrays of rings by Deblok 
et al. [38] and in a strikingly clear manner on a single mesoscopic capacitor by Gabelli 
et al. [39]. 

The second, dissipative term, is governed by a resistance R q which we call a 
charge relaxation resistance [26]: 

^ = irar (13) 

Equations (12,13) are given in the low temperature limit ksT <C hw. To understand 
the meaning of these quantities better it is instructive to consider a basis in which the 
scattering matrix is diagonal. All eigenvalues of the scattering matrix are of the form 
cxp(i£„) where ( n is the phase which a carrier accumulates from multiple scattering 
inside the cavity. Thus the density of states (4) can also be expressed as 

v = -(l/2n)J2(dCn/d(eU)). (14) 

n 

Similarly we can express the charge relaxation resistance in terms of the potential 
derivatives of phases. In the low temperature limit R q is determined by the sum of 
the squares of the dwell times divided by the square of the sum of the dwell times, 

h EJrfCnMeLQ) 2 
2e 2 (E n dCn/d(eU)y 

We now briefly discuss the charge relaxation resistance. First we note that the 
resistance unit is not the resistance quantum h/e 2 but h/2e 2 . The factor two arises 
since the cavity is coupled to one reservoir only. Thus only half the energy is dissipated 
as compared to dc-transport through a two terminal conductor. We emphasize that the 
factor two is not connected to spin: h/2e 2 results from a single spin polarized channel. 
In the single channel limit, Eq. (15) is universal and given just by h/2e 2 . This is 
astonishing since if we imagine that a barrier is inserted into the lead connecting the 
cavity to the reservoir one would expect a charge relaxation resistance that increases 
as the transparency of the barrier is lowered. Indeed, if there is a barrier with 
transmission probability T per channel in the lead connecting the cavity and the 
reservoir, then in the large channel limit, for TN 3> 1, R q is 

R q = (h/e 2 )(l/TN). (16) 



p _ "• Z^n^n/"-^ II n -s 
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Figure 3. Mcsoscopic capacitance distribution Pi,, (a) of a cavity connected via 
a single perfect lead to a reservoir. The capacitance is in units of the geometric 
capacitance, a = C^/C. The ratio rj = A(C/e 2 ) is the level spacing in units of 
the Coulomb energy. After Ref. [14]. 



In the large-channel limit, Eq. (16) is inversely proportional to the total transmission 
NT. Thus in the large channel limit Eq. (15) behaves as expected. For a one-channel 
connection to a reservoir the recent experiment by Gabelli et al. [39] indeed finds a 
resistance of 12 kil for a mesoscopic capacitor formed with the help of a quantum 
point contact in high magnetic fields in the quantum Hall regime. A brief overview of 
charge relaxation resistances in mesoscopic systems is presented in Ref. [40]. 

For a single channel connected to a cavity, the charge relaxation resistance is 
universal and given by h/2e 2 . Thus it is only the capacitance C M which needs to 
be investigated further. Gopar, Mello and one of the authors [14] and independently 
Fyodorov and Sommers [13] calculated the distribution of dwell times wp(Td) for a 
cavity coupled to a lead with a perfect one-channel quantum contact in the cases 
of time-reversal, broken time-reversal and broken spin-inversion symmetry (denoted 
by Dyson symmetry indices — 1,2,4 respectively [41]). This permits to find the 
distribution function of the capacitance (12) of such a cavity. Fig. 3 shows the 
distribution Pi^a) in the presence of time-reversal symmetry, = 1, as a function 
of a — C^/C — Td/(Td + 77), where 77 = CA/e 2 is the ratio of the level spacing and 
Coulomb energy. The distributions for larger 77 are wider, because with 77 — > 00 the 
mesoscopic fluctuations of v have a stronger effect on C M . Typically the Coulomb 
energy is much larger than the level spacing, the parameter 77 is small and the 
capacitances are close to the geometrical value (see Fig. 3). 

Next let us consider a cavity that is connected via N = 2 open quantum channels 
to the outside. The two open channels can be provided by the same contact (a single 
lead as in Fig. 2) or from different contacts as shown in Fig. 4. There are now two 
dwell times which in general differ from one another. Their joint distribution (for 
arbitrary N) is derived in Ref. [18]. As a consequence we now have a nontrivial 
distribution also of the charge relaxation resistance [27] shown in Fig. 5. For the case 
of time-reversal symmetry, 0=1, Ref. [27] finds that R q is uniformly distributed 
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Figure 4. Chaotic dot connected via two quantum point contacts with N\ and 
N2 open quantum channels. The dot is capacitively coupled to a gate at voltage 
Vn. 

between R q = h/4e 2 and R q = h/2e 2 , whereas for the broken time-reversal symmetry, 
/? = 2, the distribution exhibits a peak and reaches zero at the limiting values [27]. 
The limiting values are reached if both channels have dwell times close to each other 
(R q = h/4e 2 ) and if one of the dwell times becomes very much longer than the other 
one (R q = h/2e 2 ). The probability of such processes for (3 = 1 is higher then at (3 = 2, 
so that the distributions are stronger at R q — h/2e 2 1 h/Ae 2 . 

4. Charge noise in the presence of shot noise 

We have derived the capacitance and the charge relaxation resistance R q from 
linear response. The dissipative part of the response proportional to R q is related 
to the charge and current fluctuations of the system via the fluctuation dissipation- 
theorem. Returning to the single- lead configuration in Fig. 2, we can calculate the 
current fluctuations at contact a = 1 or at the gate a = or their correlations, 
S a p(t) = 1/2 / dt'((I a (t + t')I f j(t) + I f j(t)I a (t + t'))). Due to current conservation we 
have S(t) = Sn(t) — Son(t) = — Sio(i) = —Soi(t). For the Fourier transform of this 
noise spectrum we have S^](uj) — 2/cbT 'w 5 'C 2 R q , where we have taken the classical 
limit ftid <C ksT. Since the current is the time derivative of the charge accumulated 
on the cavity we immediately also find the spectrum of the charge fluctuations Q on 
the cavity, S e £ Q (uj) = 2k B TClR q . 

Next consider the case of a cavity with two contacts which permit exchange of 
carriers with reservoirs, as shown in Fig. 4. There N\ quantum channels in the 
QPC connect the cavity to the left contact at voltage V\ and N2 quantum channels 
connect the cavity to the right at contact 2 with voltage V2. If these two voltages 
differ, a transport state is established. At zero temperature we have a current 
/ = GV with V = V\ — V2 determined by the Landauer formula, G — (e 2 /h)T 
where T = tr6>^l2<Sli. For a chaotic cavity the ensemble averaged total transmission 
probability is T = N1N2/ '(N\ + A^). For simplicity, we consider a voltage eV > fceT 
so large that the thermal noise in this conductor can be neglected. However, the 
granularity of charge and the fact that to every incident channel there are different 
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Figure 5. Charge relaxation resistance R q , measured in units of h/e 2 , for f) ■ 
(dashed curve) and = 2 (solid). After Ref. [27] 



final channels leads to shot noise [4] with a zero-frequency noise spectrum given by 
gshot _ 2e(e 2 /h)\V\tr (li<S^li«Sli«S^l2>S). In terms of the transmission eigenvalues 
the shot noise is Sfj ot = 2e(e 2 //i)|V| Y^n^n^ - T n ). For a chaotic cavity with large 
perfect quantum point contacts N\ , N 2 S> 1 the ensemble averaged shot noise [4, 24, 42] 
is Sfj ot = 2e(e 2 /h)\V\N?N%/(N 1 + N 2 f . For single channel leads N x = N 2 = 1, the 
case of interest here, shot noise is characterized by a distribution which is given in Ref. 
[27] . The stochastic transfer of carriers through the cavity leads to fluctuations of the 
charge in the cavity as a function of time. Such charge fluctuations can build up only 
if screening in the conductor is not perfect. Therefore, consider now again a nearby 
gate which is capacitively coupled to the conductor with a geometric capacitance C. 
The fluctuating charge in the cavity induces charge fluctuations in the nearby gate 
such that the charge on the conductor and the gate is conserved. In the transport 
situation considered here, the different voltages in the leads break the symmetry 
between contacts. An important role is now played by certain selected elements of 
the full Wigncr-Smith matrix. The elements which retain information on the contacts 
from which carriers are incident are 



2m V d(eU) J 0J 



(17) 



Introducing the second quantization operators ap{E) which annihilate an electron in 
the incoming channel in lead [3, we write the charge operator [27] in the absence of 
interaction as 



ej\f = eNfo{E)al{E)a^{E). 



(18) 
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Figure 6. Distribution of the resistance R v which determines the low-frequency 
charge fluctuation spectral density of a chaotic cavity connected to reservoirs with 
two single channel leads. Dashed line is for the orthogonal ensemble, (3 = 1, solid 
line is for the unitary ensemble, (3 = 2. After Ref. [27] 



In reality the Coulomb interaction leads to potential fluctuations inside the cavity 
and these fluctuations contribute in turn to the charge. Since charge and potential are 
linearly related we are lead to describe the fluctuating potential also with an operator 
U that contributes to the net charge in proportion to the density of states v in the 
cavity. We have Qp 7 = eN^ — e 2 v\J[j 1 . On the other hand the Coulomb interaction 
dictates that Qp 7 = CUpj. Taken together, we have therefore, 

Qp-f = CUp-y — eNp-f — e 2 i>Up 7 . (19) 



Solving for \Jp 1 determines the charge fluctuations in the presence of screening, 

Q M = = e %AW (20) 



Here C p is the electrochemical capacitance, Eq. (12). Evaluation of the quantum 
statistical expectation value leads to a low-frequency charge-fluctuation spectrum 
Sqq(0) = 2e|y|C^i?„ determined by the resistance [27] 

For a chaotic cavity connected only to two single channel leads the distribution of the 
resistance R v is shown in Fig. 6. In many cases the mesoscopic distribution P(C /i ) of 
the electro-chemical capacitance C M is much narrower then that of R v , cf. Fig. 3 for 
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small rj and Fig. 6, and the fluctuations of noise Sqq are mostly due to fluctuations 



It is interesting to calculate not only the spectrum of the noise induced into 
the gate but also the spectrum at the contacts which connect the cavity to reservoirs. 
Extending the discussion of Ref. [27] Hekking and Pekola [43] find the finite frequency 
noise spectra in the quantum limit. A quasi-classical discussion of the frequency- 
dependent third cumulant of chaotic cavities is given by Nagaev, Pilgram, and one 
of the authors [44]. It turns out that in higher cumulants both the i?C-time and the 
dwell time are relevant [44] . For a related discussion of dynamical thermal effects we 
refer the reader to Reulet and Prober [45] . 

5. Relaxation and dephasing rate of a charge qubit 

To illustrate the physical significance of the charge relaxation resistance (13) and 
the resistance R v (21), we now consider a chaotic cavity coupled capacitively to a 
double quantum dot as shown in Fig. 7. The double dot represents a charge qubit: a 
single charge tunnels between the upper and lower dot. The chaotic cavity plays 
the role of a detector. With this set-up we can investigate the question [46]: if 
we use a generic mesoscopic conductor (here a chaotic cavity) how good a detector 
would this represent? The properties and suitability of special systems, like quantum 
point contacts or single electron transistors, as detectors have been widely discussed 
[30, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55] and continue to be a subject of research 
[56, 57, 58, 59, 60]. Interestingly a chaotic cavity connected to one channel leads is 
with a probability of almost one-half close to an ideal detector and if it is not ideal 
changing the shape of the cavity just a bit leads with probability of one-half to an 
almost ideal detector [46]. 

The two- level system is represented by the Hamiltonian Hdd = %b z + ^b x , 
where bi denote Pauli matrices. The energy difference between upper and lower dot 
is e and A accounts for tunneling between the dots. The full level splitting is thus 
Q = Ve 2 + A 2 . The chaotic cavity and the two-level system are coupled through 
the long range Coulomb interaction taken into account with a set of capacitances 
C\,Ci,Ci that link the charge Q\ and —Q\ on the double dot to the charge on the 
cavity. The series capacitance is C _1 = Cf 1 + C^ 1 + C^ 1 . 

If the two-level system is in a superposition state, it experiences relaxation with a 
rate T re i towards the ground state due to charge fluctuations in the cavity. In addition 
a coherent state of the two level system is dephased with a rate Td ec . 

At equilibrium the relaxation rate and decoherence rate generated by the cavity 
are determined by the charge relaxation resistance [46] 



A charge relaxation resistance exists even if the conductor permits no transmission of 
carriers. Often the interaction of a qubit and a detector is modeled by describing the 
effect of the qubit solely in terms of a modulation of the transmission amplitude of 
the conductor. Such a description predicts that a conductor with zero transmission 
probability causes no relaxation and no dephasing. But a conductor with transmission 



of R, 




(23) 
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probability zero is just a capacitor, and, as we have seen, a mesoscopic capacitor has 
a non- vanishing charge relaxation resistance. 

In the zero-temperature limit, the two-level system relaxes due to the zero point 
fluctuations of the conductor. A qubit in an excited state is a detector of vacuum 
fluctuations. If a voltage is applied to the conductor, the decoherence rate contains 
in addition a term which arises due to the charge fluctuations in the presence of shot 
noise of the detector. This additional term is proportional to the applied voltage and 
is governed by R v , 



Thus the density of states of the cavity v, its electrochemical capacitance C M , and 
the two resistances R q and R v are the properties of the mesoscopic conductor which 
determine the relaxation and decoherence rate of the two-level system. 

6. Chaotic cavity as a Quantum Detector 

The cavity in Fig. 7 can be viewed as a detector of the state of the two-level system. 
The current through the cavity depends on whether the charge Q\ is close to the cavity 
or in the dot further away from the cavity. The difference of the currents AI = h—h 
is evaluated using the Landauer formula 



where AG is the change of conductance between the two states of the double dot and 
AU = eCfj,/ ((ve 2 )(Ci — C M )) is the potential change on the mesoscopic conductor [46]. 
To measure the current, the shot noise of the cavity must be overcome. The zero- 
frequency shot noise is Su = e| V|(e 2 /27r) ^ T„(l — T n ). Therefore, a measurement 




(24) 




(25) 
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time [51, 54] r m = 45///(A/) 2 is needed for a signal to noise ratio of 1. Using weak 
coupling Ci , C 2 <C Cj one gets for the inverse measurement time [46] 

r" 1 = T m = 2^ i? m e|t/| (26) 



with the resistance 



p h (Y, n dT n /d(eU)f 

m 4 e 2(E„#n/W) 2 EJn(l-rn)' 1 ' 

Fundamentally the measurement is always slower than the decoherence, the decay of 
the off-diagonal elements of the reduced density matrix of the two-level dot. This 
implies the inequality T m < Tdec or R m < R v . This can be shown by deriving a 
Schwarz inequality for a combination of scattering matrices and the element A/21 of 
the Wigner-Smith matrix (sec Rcf. [46]) or more generally within linear response 
detector theory [55]. 

An efficient detector provides a maximum of information for a minimum of back- 
action {i.e. decoherence) on the measured system [55]. The efficiency of the detector is 
determined by the ratio r\ = T m /T dec = R m /R v < 1. The most efficient measurement 
requires that the tunneling between the two double dots is negligible, A ~ 0, and 
the temperature must be much smaller than the applied voltage k^T -C e\V\. But 
more importantly with the general formulation given above, we can now discuss the 
conditions on the scattering matrix for the detector to be ideal. Ref. [46] finds that 
the scattering matrix needs to be of block-diagonal form: channel mixing detectors 
are not ideal. Consequently the scattering matrix of an ideal detector can be divided 
into 2x2 blocks of the form 

(n) _ / -zVT^7>^+^") V?^ "-»*-> \ . . 

Each block is defined by its transmission probability T n and three scattering phases 
4 l n,4'A,n,<f>B,n- Using the definition of R v (Eq. (21) we arrive at [27] 

p h En {(dT n /dUf/(AT n (l - T n )) + T„(l - T n )[d{<f> A , n - 4> B , n )/dU) 2 ) 

{EndK/dUf 

Eq. (29) reduces to earlier results [47, 51, 52] in the infinite capacitance limit where 
C^R V in Eq. (4) can be replaced by R v {ve 2 ) 2 . In Eq. (29) the derivatives dT n /d(eU) 
determine the sensitivity of the cavity to a potential variation AU. A high sensitivity 
implies a fast detector (a short measurement time). However, a high sensitivity also 
implies a fast decoherence rate proportional to R v . From this we see that decoherence 
and measurement speed are closely related quantities. 

Demanding that R m and R v are equal, requires d(j>A,n/dU — d<pB,n/dU — 0. As 
pointed out by Clerk, Girvin and Stone [55] , this condition is not connected to any 
physical symmetry of the system, but can be understood as a condition that the 
detector should not transfer information into the phases of the scattered electrons 
since these are not measured. However, we can require that each individual phase- 
dependent term vanishes separately. d<j> b ,n I 'd(eU) = in Eq. (29) can be achieved 
by requiring that the scattering Hamiltonian must obey time-reversal symmetry. The 
phases d(j>^ n /d(eU) — vanish for detectors that obey a spatial inversion symmetry 
V(x,y, z) = V(x, y, — z). (Here we assume that conduction is along z and x,y are 
transverse coordinates). Interestingly Ref. [46] finds that in the multichannel case 
N > 1 another condition is needed! The equality R m = R v gives us a condition which 
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Figure 8. Efficiency distribution of an ensemble of chaotic quantum cavity 
detectors: orthogonal ensemble (top panel), unitary ensemble (lower panel) for 
single channel (N\ = N2 = 1) and double channel (Ni = N2 = 2) point contacts. 
After [46] 



is of statistical origin. The total conductance of the detector is a sum of one channel 
conductances that have independent uncertainties. The statistical uncertainty of this 
sum is minimized if [46] 

5y£?5y - C(U»U). (30) 

with a function C(U) > that does not depend on the index n. In the WKB limit 
we have d/dE = —d/d(eU) and Eq. (30) can be interpreted as differential equations 
for the transmission probabilities T n (E). The solutions are all of the form T n = 
(1 + e -(F(E)-F(E n ))yi with dF / dE = c (The function F is therefore monotonously 
increasing). The only difference allowed between the different probabilities T n is the 
offset energy E n . Transmission probabilities of the type (30) occur automatically if 
the scattering problem is separable due to a potential of shape 

V(x,y,z) = Z(z) + W(x,y). (31) 

This occurs in particular for the case F = 2ttE/lj z with a symmetric harmonic 
scattering potential Z{z) = Vq — muilz 2 /2 which is the saddle point potential of 
an ideal quantum point contact. 

To illustrate the role of the condition Eq. (31) we now consider chaotic cavities as 
detectors. The condition Eq. (31) states that a geometry with a separable potential 
V(x, y, z) — Z{z) + Y(x, y) is favorable to obtain an efficient detector in the case of 
more than one open channel. It is clear that a chaotic cavity violates this condition. 
The efficiency rj — r m /r„ = R m /R v of a chaotic detector with two open channels in 
each contact is expected to be much smaller than a chaotic cavity connected to single 
channel contacts for which Eq. (31) plays no role. 

Using the distribution of the elements of the Wigner-Smith matrix (the Lagucrrc 
ensemble [17]) we find the probability distribution of the measurement efficiency 
V = Rm/Rv in the orthogonal (time-reversal symmetry) and unitary ensemble (broken 
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time- reversal symmetry) shown in Fig. 8. The distribution for Ni = N2 = 1 in 
the unitary ensemble can also be calculated analytically, P{rj) = l/y/rj(l — rj). The 
distributions for the other cases were obtained by numerical integration. Interestingly 
Fig. 8 shows that despite the absence of inversion symmetry a chaotic dot with open 
single channel contacts is with high probability an efficient detector! It is clearly 
visible that chaos reduces strongly the efficiency of the measurement device as soon as 
more than one channel contributes to the electric transport. Compared ti the time- 
reversal symmetric case (solid curve in the upper panel) the reduction due to a broken 
time-reversal symmetry (solid curve in the bottom panel) is much less pronounced. 



7. Quantum pumping 

Quantum pumping exploits the sensitivity of quantum interference to variations of 
parameters of the scattering matrix. The variation of two parameters X\^{t) — 
Xi^ cos(wf + 0i, 2) oscillating with the same frequency but out of phase generates 
a dc-current [61, 62, 63, 64, 65, 66, 67, 68], if the oscillating scatterer is in a zero- 
impedance external circuit, or a dc- voltage [69] (if the pump operates in an external 
circuit with finite impedance). Oscillating parameters X12 are due to, e.g. voltages 
applied to the gates that form the shape of the dot. Here we briefly discuss quantum 
pumping, see Fig. 9, when the dc voltage generated by the pump is measured, similar 
to the experiment by Switkes et al. [69] . For slow variations of the shape of the dot 
we can use the scattering matrix S(e) and later include Coulomb interactions self- 
consistently[62, 70]. The Coulomb potential is found from the condition of charge 
conservation. Using a global gauge transformation it can be shifted into the phases 
of the the scattering states. However, such a transformation in itself does not pump 
electrons. The potential can not be considered as an additional independent pump 
parameter similar to X(t). 

We consider a two-terminal quantum dot with two weak potentials X\^ varied 
with the same frequency lu. In order to derive a formula for the dc voltage V, we use 
a simple model for the quantum dot and the two electron reservoirs, see Fig. 9 a. The 
dot and the reservoirs 1, 2 are connected to a screening gate via capacitances C and 
Ci.2- Following Rcfs. [71, 72], we introduce the emissivity edq(a)/dXj and edq(a)/dE, 
which is the charge that exits the dot through a point contact a (a = 1,2) when the 
parameter X or chemical potential E is changed: 

- J-trS+l — - — trSti ^ f 32 ) 

dx ~2™ tr6lQ ax' be ^2m tTi,la dE- [62) 

The total current I\ flowing through contact 1 is 

i— 1 

where g = tr5^1nSl2 is the conductance of the quantum dot. A similar expression 
determines 12- For slow variations of X we can restrict ourselves to first order time- 
derivatives only. Since Ij = CjdVj/dt, in the limit w <C G/C\ t 2 we have I\{t) — 77Z2 (i) , 
where 77 = C1/C2 is a numerical coefficient describing the capacitive division between 
the two reservoirs. The dc voltage for small periodic parameters with amplitudes Xi^ 
and phase difference <ft = <Pi — 4>2 is V — (/iw/4e)XiX 2 £sin0, where [73] 

g /l fl(g(l)-T?g(2)) \ _ _d_ fl d{g{l)- m {2)) Y 
8X 2 \g dX 1 ) dX 1 \g 8X 2 )_ 



1 

v = - — 



(34) 
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Figure 9. (a) Schematics of the voltage measurement set-up: Quantum dot 
and reservoirs with potentials Vi(2) are capacitively coupled to the ground via 
C and Ci(2); (b) Weak time-dependent potentials -Xi,2(t) generate pumped dc 
current (voltage); (c) The pumped voltage V is proportional to the area of the 
contour in parameter space of the voltages Xi ; 2(t) = -Xi,2 cos(uit + ^1,2) (shown 
for <fii — <f>2 = vr/2). 



The derivatives with respect to X\ and X 2 in Eq. (34) are taken at constant values 
of the chemical potential E of the reservoirs. However, when the Coulomb interaction 
is accounted for, it is the sum of chemical potential and local self-consistent Hartree 
potential U, which is constant. Therefore, in the presence of Coulomb interactions the 
expressions for the transmitted charge should take this into account [70] . In the limit of 
weak intra-dot interactions, e 2 /C <C A, one still can apply the non-interacting theory 
and take all derivatives with respect to shape- varying potentials X\^ at fixed value of 
E. On the other hand, in the strong electron-electron interaction limit, e 2 /C 3> A, 
which is usually relevant in experiments, transport occurs on the background of almost 
constant charge q in the dot. Following Brouwer [62] we can relate derivatives with 
respect to X at constant E and q and in this way include a Hartree potential of 
arbitrary strength. 

First we note that the applied potentials X, the chemical potential E and the 
charge of the dot Q are not independent variables. If we denote the derivative of some 
function F with respect to X at constant Y as (dF/dX)y, then 

which allows us to write 

d\ ( d \ (dE\ ( d \ {dg/dx) E a 

The derivatives in the last equation in Eq. (36) are found as follows: the denominator 
is evaluated in the absence of pumping (constant X) as a balance of particle current 
and displacement current from the dot, (%2 a p G a fj{u)) +iu;C)dE = iuuedq^, so that the 
derivative is dq u /dE = {C/C^trN. The numerator is evaluated at constant E, so 
the dependence of S on the pumping parameter X is crucial. Unlike the dc-pumped 
current, found by Brouwer to be quadratic in pumped strengths [62], the ac-component 
I u is linear in X and can be found by expansion of S([X(t)]) in X to first order. As 
a result, dq^/dX^ = tvMx- If analogously to Eqs. (1,2) we introduce the matrix 
M x - {l/2iri)S^d x S then 

8E_ C^trfSx . . 

OX C tvN ' [ ' 




Figure 10. Mcsoscopic distribution P(v) of pumped voltage v for Ni = N% = 1 
with and without time-reversal symmetry (solid and dashed curves respectively). 
The limits of strong, C/e 2 <C A (main figure, after Ref. [73]), and weak 
interaction, C/e 2 ^> A (inset), are similar to each other. 



In the non-interacting limit C — > 00 only the first term in the r.h.s. of Eq. (36) 
survives. In the opposite limit C — > of a realistic quantum dot the charging energy 
is large, e 2 /C » A. Then one finds that the charge of an open dot essentially remains 
constant during the pumping cycle, I\(t) = —hit) for all time t. As a consequence, 
the pumped voltage v loses its dependence on the capacitive division 77, the ratio 
C^/C — * 1 and therefore the derivatives with respect to X are 

(dx) q = (dx)E- t -^d E . (38) 

Using Eqs. (36,37) and using the emissivity dq/dE (32) we reformulate Eq. (34) for 
arbitrary interaction. The distribution of Mx in the numerator of Eq. (38) is known 
[18], so that one can perform numerical integration for few-channel dots or use the 
diagrammatic technique [74] to find the distribution of relevant quantities, e.g. the 
pumped voltage v. 

The pumped voltage v is zero when averaged over an ensemble of dots, and the 
width of the distribution P(v) is expected to diminish with growing N, similarly to the 
distribution P(i) of the pumped current [62]. Ref. [73] finds that in the multi-channel 
limit the distribution has a Gaussian shape with r.m.s. (u 2 ) 1 / 2 = (2/ttN 2 ), since 
the fluctuations of conductance become small and the conductance takes its classical 
value, g — N/2. However, in the few-channel limit Ni = N 2 = 1 the mesoscopic 
distribution is wide, see Fig. 10. The distributions in the limits of weak and strong 
Coulomb interaction, treated with the help of Eq. (38), are similar to each other (in 
the weak interaction limit we take symmetric capacitive division, 77 = 1). Particularly, 
for (3 = 2 numerics gives almost identical curves. Thus the self-consistent internal 
potential does not change distributions significantly, a conclusion that, for pumped 
currents, was reached previously in Ref. [62] . 



8. Shot noise of photon excited electron-hole pairs 

In this section we consider the charge fluctuations and the zero-frequency current noise 
generated by oscillating potentials applied to the contacts and the gate of a multi- 
channel chaotic sample with M > 4 leads (see Fig. 11). Several contacts are subject 
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Figure 11. Schematics of a noise measuring set-up. A four-terminal (M = 4) 
chaotic dot is subject to oscillating potentials V a (t) at contacts a = 1,2 and 
coupled to a gate with a time-dependent potential Vb(t), via a capacitance C. 
The internal potential of the dot is U(t). Tunneling elements (shaded areas) 
partially close the dot. 



to periodic ac-potentials V a (t) = V a cos(uit + (j) a ) at the same frequency u> but possibly 
with different phases (j) a . We consider the case eV a <C fiu, where the amplitude of 
the applied potentials eV a is small compared to the modulation quantum h\o. The 
oscillating potentials excite the electron gas in the contacts by creating electron-hole 
pairs. These photon-excited carriers are then reflected at the sample or transmitted 
through the sample into another contact. There is no dc-current linear in voltage. 

We find effect of Coulomb interactions on the noise through such a multi-channel, 
ac-biased (partially) open quantum dot for arbitrary frequencies u) and interaction 
strength, see Fig. 11. Since in an experiment with dots or similar structures [75, 76] 
frequencies u> might be comparable to the inverse dwell time, one can not restrict 
oneself to first energy derivatives only. Consequently, the Wigner-Smith matrix is not 
as useful for arbitrary u>, since we have to consider correlations of electrons at quite 
different energies. However, in the limit N ^> 1 for open, or high total transmission 
trf 3> 1 for partially open cavities we can investigate correlations at arbitrary u>. 

As emphasized in the Introduction, the Coulomb interaction is very important 
for frequency-dependent transport. In the non- interacting theory the noise of photon- 
assisted carriers was discussed by Lesovik and Levitov[77] and Pedersen and Buttiker 
[78] for energy- independent scattering matrices. This is a good approximation if the 
frequency u> is sufficiently small, uot^ <C 1. The sample-specific results of Refs. [77, 78] 
are in a good agreement with an experiment by Reydellet et al. [75] in quantum 
point contacts. However, if an experiment is performed at high frequencies, such that 

~ 1 , the non- interacting theory provides a result which is not gauge- invariant [79] . 
In reality, the frequency-dependence should be treated by taking Coulomb interactions 
into account. In the set-up shown in Fig. 11 we consider the current correlations in 
the leads A, \i using scattering theory for chaotic systems [2, 4]. 

The rate of photon-generated pairs is small [80] . If the pairs are split into different 
contacts, they generate current noise. For a single ac-voltage this noise corresponds 
to electron-hole correlations of elements of a pair. For several ac-bias voltages an 
adjustment of the phase shift between them allows one to maximize (minimize) noise 
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and extract information about the scattering matrix [80] . This is why it is interesting 
to find current noise correlations Sx^ between A-th and fi-th lead for a multi-lead 
(partially) open quantum dot. Similar results were obtained by Samuelsson and the 
authors in Ref. [79] for completely open dots. 

Below we first discuss shot noise for open quantum dots and show why it is 
necessary to take the Coulomb interaction into account. Subsequently we discuss the 
noise correlations to leading order in trf for partially open dots, where the total 
transparency is still large trT 3> 1. In this limit the mesoscopic fluctuations of the 
noise S can be neglected, so that the result gained from mesoscopic averaging is 
representative. Corrections to the noise due to the symmetry appear in sub-leading 
orders, so we do not distinguish here between (3 = 1,2. The coherent open dot with 
N channels is fully characterized by its N x N scattering matrix S. Scattering is 
spin-independent and the results given below are presented for a single spin direction. 
The dot is assumed fully coherent and the effects of inelastic scattering and dephasing 
are discussed in Ref. [79]. We use units e = h = fcs = 1. 

We start by expressing the sample-specific noise Sx^ in terms of energy-dependent 
scattering matrices 5(e), amplitudes of applied voltages V a , their phases 4> a and the 
electron distributions in the leads f a (s)- Ref. [80] finds, 

2 

Sx„ = |^Re ]T tr (A a0 (\,mu)A 0a (»,O)) (s m , V* - V a V p 6 mil S*»-*^) , (39) 

m. a/3 

where A(A, e) = 1\ - 5 t (e)l A 5(e). In the non- interacting low-frequency limit, such 
that S is energy- independent, the ensemble averaged noise is [79] 

S °^ = ^( Nt * V 2 -\trVe^), (40) 

where we introduced the averaged conductance of the open quantum dot 
G\fj, — (N\5\p — NxN^/N) and diagonal matrices of the amplitudes V = 
diag(Vili, VmIm) and phase shifts = diag(0ili, 0m1m)- This result is gauge- 
invariant, that is an arbitrary uniform shift of all Fourier components of potentials 
Vclu = V a exp(i<j) a ) by some (complex) potential Uq does not have any physical 
consequence. We notice here that although formally the 2nd term in Eq. (40) is 
of order 0(1/N 2 ), the traces in the numerator can be of order N, so that both terms 
may be of the same order. For example, when a global shift of potentials is applied 
both terms are equally important. 

However, if we are at sufficiently high frequencies, uTd ~ 1, the energy-dependence 
of S is important, and the noise spectra are given by 



2G Am ItrFe^l 2 
'a m - + N 2 ul + ( WTd )-2 



These spectra are not gauge-invariant because of the 2nd term. This is a consequence 
of the fact that we neglected the internal potential U of the dot and the effect of 
external gates with potential V cos(u;i + 0o) capacitively coupled to the dot. Therefore 
now we consider explicitly the effect of the induced internal potential U u . The sample- 
specific formula (39) for the non-interacting system can be used again if we shift all 
potentials V u — > V u — U u , including unbiased leads. 

To proceed we find the mesoscopically averaged potential U u for a partially open 
dot. The mesoscopic fluctuations of the potential (9) around the average are small 
due to the large parameter trf ^> 1, so we are justified to use if we arc interested 
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in current correlations only to leading order. As a consequence, any effects related to 
time- reversal symmetry are neglected both in U u and, in the end, in the noise S Xfl . 
At high frequencies, wr d <~ 1, we have to fully account for the frequency dependence 
of U u . The leading order of Eq. (A. 4) for the pair correlator of energy-dependent 
scattering matrices [5, 81] allows us to hnd the average of the potential U u and its 
derivatives, 

U "- V ° e + 1 - icoRC, { trf V ° ) ' 8V a , u 1 - icoRC, trf ' ( ' 

where we introduced R = 1/trT, which corresponds to the ensemble-averaged charge 
relaxation resistance (R q ) given in Eq. (13) for a partially open dot, see also Eq. 
(16). For an open dot, T = 1 the result corresponds to that of [79]. We see that the 
dwell time Td = l/(trTA) of matrix correlators is replaced by the i?C-time of charge 
relaxation inside the dot, r = (C^R. This substitution usually occurs in leading order 
in N. The ui — > 0, T = 1 limit of Eq. (42) could be easily obtained from Eq. (10) using 
the Wigner-Smith matrix and ensemble averaging (using results of Ref. [18]). From 
Eq. (42) for r = 1 and a shift of all potentials V u<a as described above, we obtain 
with a little algebra [79] : 

The second term in this equation is obviously gauge invariant and should be contrasted 
with that in Eq. (41). Experimentally, one usually has ujt <C 1 because the Coulomb 
interaction is sufficiently strong, and the 2nd term in Eq. (43) vanishes. Thus the non- 
interacting low-frequency limit (40) is recovered. The fact that the limits of strong 
interaction at lot <C 1 and weak interaction at uot^ <C 1 provide the same result 
explains why the results of the non-interacting theories of Rcfs. [77, 78] are in good 
agreement with experiment [75]. At sufficiently high frequencies, which are in the 
range of modern experiments, the difference will become apparent. 

If one now considers a partially open quantum dot with channel transparencies 
r.; 7^ 1, a similar treatment is possible. For simplicity we consider the case when 
all channels have the same transparency Ti = T. One has to use Eq. (A. 6) for 
the correlators of scattering matrices. As above, the external gates with potential 
Vo(t) = Vq cos(uit + 4>q ) are capacitively coupled to the dot via a capacitance C. To 
demonstrate the gauge-invariance of the final result for the noise, we assume that the 
current-measuring leads A, fj, are also biased. Then the non-interacting noise is given 
by an expression analogous to Eq. (40): 

2r (,r x ™ x na _ 2 , T ^ l2 , , r(i-r) 



*v(r) = Jr. U §Xfi E^M {NtI V 2 _ |tr Vei , f) + . 



J K 1 1 ; N 2 LU 

N x S x ,-^p^j \t T Ve i *-NV x e i <t»\ 2 +...\ x ^, (44) 

where the last term (...) is obtained from the 2nd by swapping A <-» fx. If we now 
take into account the frequency dependence of the scattering matrices similarly to the 
completely open dot, and use Eq. (A. 6) we find: 

^ - jfc k - ^) v ' - '* v ^ + it, T + (jy 



, rq-r) / 2n x n, 



tr Ve* - iurJVV e* i 
NV x e 



1 — iLOT 



2 



(45) 
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As an application, consider the cross-correlations in the unbiased leads, ^ / /i and 
V X = V I1 = 0, when V = 0: 

Notice that the interactions can not change the negative sign of cross-correlations. 
For sufficiently transparent barriers, T > 2/3, the Coulomb interactions diminish 
the absolute value of the cross-correlation (46). This is easy to see from the 
dependence of the RC-time on the geometrical capacitance, t _1 = trf(A + e 2 /C). 
Indeed, the interactions in the form of the Hartree potential lead to additional 
(displacement) currents, so that their fluctuations "damp" particle current fluctuations 
[79]. On the contrary, for T < 2/3 the 2nd term in Eq. (46) becomes positive and 
Coulomb interactions enhance the noise. Indeed, for T — > the response of the dot 
becomes capacitive rather then resistive (see the denominator in Eq. (42)), and the 
displacement currents fluctuate out of phase. 

Next, similarly to van Langen and Buttikcr [82], we consider the exchange 
correlation P ex (0). This is the difference of the noise S34 with bias Vi jW = V 2 *w = 
Ve^l 2 and of the sum of the noises for Vi = 0, V2 = V and V\ = V, V 2 = 0. For a 
four-channel dot we find 

P cx (0) = S M {Ve^/\Ve-^) - S 34 (V,0) S 34 (0, V) = (47 ) 

lOW 1 + U> Z T Z 

For a dc-biased dot, the inversion of the sign of the exchange correlation at T = 2/3 
was found previously [82]. Notice that Coulomb interactions enhance the absolute 
value of the exchange correlation (47) . 

9. Discussion of related works 

In this section we briefly allude to alternative approaches to treat the Coulomb 
interaction and discuss the applicability of the results presented above. In this paper 
we used the RPA approach to the Coulomb interaction which is built into the scattering 
matrix approach by Buttiker and co-workers [71, 72, 83, 84]. Another, often favoured 
path is the Hamiltonian approach. For non-interacting electrons it is completely 
equivalent to the scattering approach [85, 86] and which method one uses is a matter of 
taste. Aleiner, Brouwer, and Glazman [5] demonstrated how the Coulomb interactions 
can be considered in the Hamiltonian approach. The interaction Hamiltonian has a 
hierarchical structure: its universal part depends only on the number of electrons 
inside and the energy -scale e 2 /C The non-universal terms due to non-uniformity 
of the potential in the dot [5, 87] are small, <C A. Writing the Hamiltonian of the 
total system as a sum of Hamiltonians of the closed interacting dot, of the leads, and 
the coupling term, one can express transport quantities and proceed with a Green 
function formulation [5, 88]. Then the Coulomb potential is usually considered in 
several theoretical limits: cither for arbitrary number of channels N, if the interaction 
strength is small [89, 90, 91, 92], or in a diagrammatic expansion in 1/N <C 1 for 
arbitrary interaction strength [93, 94, 95, 96]. One can see that the total transmission 
tr r is an essential parameter of the problem. 

It is known from the literature that the mean-field treatment, which renormalizes 
the interaction line, can not lead to such effects as, e.g. charge quantization in 
blockaded quantum dots [21]. However, in the multi-channel limit N » 1 the exchange 
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terms can be neglected, and the Hartree potential of the dot provides the leading 
effect of the Coulomb interaction. Therefore, one expects that in a crossover from 
multi-channel systems to the Coulomb blockade regime, the exchange contributions 
may become important. Only very recently did Brouwer, Lamacraft and Flensberg 
[88] consider a unified approach using a Keldysh technique and scattering matrix 
theory. They demonstrated that the exchange diagrams may become important when 
the number of conducting channels N becomes small, and so for the dot with poor 
transmission a Hartree treatment is not enough. For the variance of the electro- 
chemical capacitance C M and the density of states v(e), as well as for the variance of 
the charge pumped by a quantum pump, the corrections due to exchange were found 
to be 0(1/N) for N ^> 1. The self-consistent theory by Biittiker and co-workers is 
reproduced for open coherent dots to leading order in 1/N. Inclusion of exchange 
terms to arbitrary order is a challenging task, but the authors [88] argue that the 
expansion in 1/N is still useful even for two-channel systems (N = 4 if spin degeneracy 
is taken into account) and the proposed diagrammatic technique in 1/N can still be 
used for few-channel dots. As another example of an application of a 1/iV-expansion 
to few-channel open dots, we mention recent work by Vavilov, DiCarlo and Marcus 
[97]. There the strong Coulomb interaction was treated as a self-averaging Hartree 
potential and was transformed into the phases of scattering states in the leads. Good 
agreement with experiment is reported for the variance of the photovoltaic current, if 
the leading in N term is used when N = 2. 

To summarize, we expect that the self-consistent treatment of the Coulomb 
interaction within the scattering approach is still qualitatively correct down to 
N = 1 — 2 channels. Measurement of transport involving sufficiently high frequency 
oj would clarify the question of the importance of the exchange terms [39, 98]. 

10. Conclusions 

In this work we have focused on transport problems in which dynamic charge 
fluctuations play an essential role. Chaotic quantum dots provide an example of 
a generic conductor which serves to illustrate the basic role which the long range 
Coulomb interaction plays in problems of this type. The isotropic chaotic scattering 
allows us in many cases to treat interaction effects with just a single potential. As a 
consequence we can simplify the theoretical discussion to an extent which would be 
unrealistic and inappropriate for less generic conductors. We have pointed out that at 
low frequencies the charge fluctuations can be described with the help of a generalized 
Wigner-Smith matrix in which we replace the energy derivative with a derivative with 
respect to the local electric potential. Together with a self-consistent treatment of the 
charge response this leads to a charge operator, which determines capacitances, kinetic 
inductances, the weakly non-linear current voltage characteristics, charge relaxation 
resistances and the low frequency noise power of charge fluctuations. We have shown 
how this approach is applied to specific examples, like the rate of relaxation and 
decoherence of a charge qubit near a chaotic cavity, or quantum pumping. We have 
illustrated how the approach is applied to transport problems in which the frequency is 
not small, discussing the shot noise of photon-excited electron- hole pairs. The unified 
description of these problems demonstrates the generality of the approach, and makes 
it a useful tool for future investigations. On the experimental side there is an increasing 
interest to explore dynamic phase-coherent transport. This interest derives not only 
from the fact that the dynamics of small structures is largely an unexplored area, but 
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also from the desire to use quantum coherent structures to perform useful tasks, like 
quantum information processing, and to perform these tasks as fast as possible. 
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Appendix A. Matrix correlators for energy-dependent matrices with 
non-ideal leads 

In this Appendix we present several correlators of the scattering matrix S of chaotic 
partially open quantum dot, when the contacts have transmission 1^ < 1, i = 1, . . . , N. 
These correlators are used to find conductance auto-correlations and noise through a 
biased dot. The imperfect transmission of the contacts is characterized by an N x TV 
reflection matrix S for which we take the simple form 

S = (S) = (1-f) 1 / 2 , (A.l) 

where T is an N x N diagonal matrix containing the transmission coefficients Tj € [0,1]. 
In order to describe coherent energy-dependent scattering with nonideal leads, we 
use the following method. An N x N matrix U(e) of a fully open quantum dot is 
found using stub model [18, 99, 100]. Energy-independent reflection from the contacts 
is modeled by S, but since the carrier can enter the dot with probability T 1 / 2 , its 
scattering is described by [86, 101] 

S(e) =S- SS(e) =S- T^ 2 U(e) (l - (1 - r) 1/2 W(e)) ~* T 1 ' 2 . (A.2) 

The matrix correlators below are presented for the case of fully broken time-reversal 
symmetry, (3 — 2. When this symmetry is present, one should add terms which 
correspond to permutations of indices, so that e.g. for 2-matrix correlator, a Cooperon 
contribution is added, and for the irreducible four-matrix correlator (Hikami box) each 
term containing four (5-functions should be supplemented by 7 other terms, obtained 
by permutations of 3 pairs of indices. 

The basic element of diagrammatic technique is a 2-matrix correlator, an 
ensemble-averaged product of two scattering matrix elements [81]. For convenience we 
introduce a shorthand notation D(e, e') = l/(tr V— i(e— e')), where e, e' are normalized 
by 2-k/A. The correlator (SSSS*), shown on Fig. Al, is 

((SS) t3 (e)(SSy kl (e')) = r<riD(e,e') (fofy + ^(1 - T t )(l - Ti)D(e, e')Mki) ■ ( A - 3 ) 

Notice that the 2nd term is actually sub-leading in tr f, but we keep it for the reasons 
to be explained later. Using Eq. (A. 3) we find 

(S l3 (e)S* kl (e')) = -r0(l-r,)(Sy4z +T i TiD(e,e')6 ik S j i 

+ v/(i - r,)(i - r ; )r 4 r^ 2 ( e , e ')^fe- (A.4) 

The leading order term ( first line) in Eq. (A.4) was previously found in Refs. [5, 81]. 
The expressions for correlators of n matrix elements for n > 2 arc more complicated, 
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Figure Al. Top: Ladder diagrams for leading order T>o contribution to the full 
2-matrix correlator V = ((5S)(5S* )). Bottom: Leading and sub-leading terms 
for "D. The 2nd term contributes, e.g., to var G. 



and below we present a 3-matrix correlator of SSSS*SS* to leading order: 
((5S) ij (s)(6S)t lh (e[)(6S)l 2l2 (e 2 )) = r^D^e'JiJ^.^) 

X (^jV 1 _ r j^2^fe2h^fclj + r iV /l - f^fc^fc^j) . (A.5) 

The 4-matrix correlator is equal to the reducible (disconnected) part and the Hikami 
box denotedB for the irreducible 4-th moment: 

((55) <1J - 1 ( e i)(55) <aJ - a ( e2 )(<^)|E lIl ( e i)(55)^ Ia (4)> = " + r <1 r ia r Il r Ia 

x (D(ei,s' 1 )D(e2,e' 2 )S ilkl 5 jl i 1 Si 2 k 2 Sj 2 i 2 + D(ei,e' 2 )D(e 2 ,s' 1 )6i 1 k 2 Sj 1 i 2 Si 2 k 1 S j2 i 1 ) , (A. 6) 

■ = r il r i2 r il r i2 D( ei , £ ;)D( e2 ,4)£)( £l ,4)£>( £2 ,4){^ lfel 5 jli2 (5 i2fe2 ^ 2;i 

x ([tr (r 2 - 2r) + i( £l - 4 + e 2 - 4)] + (1 - rjJD-^ei.ei) 

+ (i - r^ir 1 ^ 4) + (i - rjjD-^ea, 4) + (i - r^zr 1 ^^)) 



+ Si^Sj^s^Shtt ([tr (r 2 - 2r) + i( £l - 4 + £2 - 4)] + (1 - ryD- 1 ^^;) 

+ (1 - T^D-^eue'J + (1 - riJIJ-^ea.ei) + (1 - T l2 )D- l {e 2 ,e' 2 )) 

+ (ShkJhkzShjJnhO- - r)- 2 /2 (i - r); 1 / 2 + ^ 1J2 (5 l2 fc 2 ^ 2 fc 1 (5 (lJ1 (i - r)V 2 (i - r), 1 / 2 ) 

x [ j d- i ( £i ,4)+^- i ( £2 ,4)] + (w^w^a-n^a-H, 1 / 2 

+ 5 ilj2 Wh/b^ 1 (l-r)^ 2 (l -r), 1 / 2 ) [D-Hei.e'aJ + ^^.ei)]}. (A.7) 

We point out that the 4-matrix correlator does not simply reduce N to tr T or 
trT 2 , when the dot becomes partially open. As a result, Hikami box found in 
Appendix A of Ref. [81] for partially open dots reproduces Eq. (A.7) only for 
r = 1. One can use the correlators (A.4-A.7) to find, e.g. the two-terminal 
dimensionless conductance G through a quantum dot and its correlation function 
((G(e)G(e')}} = (G(e)G(e')} - (G(e))(G(e')> at different energies. To do this we 
express conductance in terms of a traceless matrix A as 

G 12 (e) = tr (St( e )i l5(e )i 2 ) = ^1 _ tr (&{e)AS( S )A), A = ^ - ^1 2 . (A.8) 
Averaging, using 2-matrix correlators, yields the energy-independent conductance: 



tr 2 fA trf 2 A 2 
(G(e)) = trTA 2 - — - - r-g-fy.i 



trT 



trT 



(A.9) 
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Figure A2. Diagrams for the evaluation of ((G(e)G(e')}): (left) A single diffuson 
T> connects one pair of matrices, the others being averaged to (<S(<S*)); (right) V 2 
gives diagrams of the same order. For (3 = 1 an additional contribution comes 
from the substitution V — > C in these diagrams. 



To find {{G(e)G(e'))) we have to calculate 2 diagrams shown in Fig. A2. 

((COW))) = f f 2 Tf«rA'(l-f)f« + tt -A-(l-f)f y'Af V A1 

U W V ;// /3 \ tr 2 r+(e-e') 2 tr 2 r + (e - e') 2 / 

Notice that for the left diagram on Fig. A2 we have to use a sub-leading term from 
Eq. (A.3), which gives the tr 2 A 2 (l-f)f term in Eq. (A.10). The doubling of the left 
diagrams comes from possibility to have correlators on the top and bottom. Due to the 
convenient representation of G in terms of traceless matrix A, the irreducible averages 
(A. 5) and (A. 7) do not contribute to the correlations. At f = Y\,Ni = N 2 — N/2 
and £ = s' we reproduce the result var G — (2 — 2r + T 2 )/ (8/3) obtained by Brouwer 
and Beenakker [74], and at T = l,s ^ e' we have var G = (l/8/?)/(l + (e - e') 2 /N 2 ) 
[81, 102, 103]. 

To obtain Eqs. (43-45) we need a correlator (A. 6) combined in the form of 
A a| g(A, e)Ap a ((i, e'), the matrices A are defined after Eq. (39). One can also use Eq. 
(A. 6) in the limit oj — > to obtain the noise in a dc-biased dot [82]. 
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